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Abstract 

In this paper, we have discussed the gravitational collapse and 
expansion of charged anisotropic cylindrically symmetric gravitating 
source. To this end, the generating solutions of Einstein-Maxwell field 
equations for the given source and geometry have been evaluated. 

We found the auxiliary solution of the filed equations, this solution 
involves a single function which generates two kinds of anisotropic so¬ 
lutions. Every solution can be expressed in terms of arbitrary function 
of time that has been chosen arbitrarily to fit the various astrophysi- 
cal time profiles. The existing solutions predict gravitational collapse 
and expansion depending on the choice of initial data. Instead of 
base to base collapse, in the present case, wall to wall collapse of the 
cylindrical source has been investigated. We have found that the elec¬ 
tromagnetic field is responsible for the enhancement of anisotropy in 
collapsing system. 

Keywords: Cylindrical Symmetry; Gravitational Collapse; Electromagnetic 

Field. 
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1 Introduction 

The study of gravitational collapse started from the pioneer work of Chan¬ 
drasekhar (1936). Later on, Oppenheimer and Snyder (19390investigated 
that the spherically symmetric homogenous dust collapse leads to the forma¬ 
tion of black hole (BH). Initially, it was argued that the homogeneity and 
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spherical symmetry of the collapsing model are responsible for the formation 
of BH. However, it was found (Lemaitre 1933) that inhomogeneous dust col¬ 
lapse would end as shell crossing and shell focusing singularities. After many 
years, it was proved (Joshi 1997) that the shell crossing singularity is naked, 
while shell focusing singularity might represent BH, depending on the choice 
of initial data of collapse. Hence, it was concluded that the homogeneity of 
the collapsing model is not a sufficient condition for the formation of BH. In 
order to generalize the collapsing matter, it becomes necessary to study the 
collapse of more realistic matter with non-vanishing pressure. 

Misner and Sharp (1964) studied perfect fluid collapse and found BH 
as the end state of gravitational collapse. Also, Herrera and Santos (1997) 
explored the properties of anisotropic self-gravitating spheres and discussed 
their stability using the perturbation method. Herrera and his collaborators 
(Herrera et al.2008a, Herrera et al. 2008b, Herrera at al. 1989, Herrera et 
ah 2009, Herrera et ah 2010, Herrera et ah 2012) have discussed the sta¬ 
bility and applications of anisotropic solutions to stellar collapse. Recently, 
Glass (2013) has formulated a generating solution of anisotropic spherically 
symmetric solutions which reveal either expansion and collapse depending 
on the choice of time profile of the solutions. This work has been extended 
for plane symmetric anisotropic source and charged anisotropic sphere by 
Abbas (2014a, 2014b). The present paper is the cylindrical version of the 
these papers. 

The current observational evidences of gravitational waves (through the 
detectors LIGO (Abramovici et ah 1992) and GEO (Luck and GEO 600 
Team 1997) have increased the interest to study the gravitational collapse in 
cylindrically symmetric systems. The spherical systems are simple and do not 
provide non-trivial examples in the generic gravitational collapse. Therefore, 
the study of cylindrical collapse is much important as compared to spherical 
systems. Some numerical studies (Piran 1978) provide the generation of grav¬ 
itational waves from cylindrical collapse. These results have been extended 
analytically (Nakao and Morisawa 2004) to study gravitational waves during 
cylindrical gravitational collapse. Sharif and Ahmad (2007) generalized this 
work for two perfect fluid cylindrical collapse and discussed the generation of 
gravitational waves. Di Prisco et ah (2009) studied the shearfree gravitational 
collapse of the anisotropic fluid in the cylindrically symmetric spacetime. 

The applications of the electromagnetic field in astronomy and astro¬ 
physics is an active research domain. A lot of work has been devoted to 
discuss the collective effects of electromagnetic and gravitational fields. Till 
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now, there is a little progress about the effects of electromagnetic field on 
gravitational collapse of stars. Thorne (1965) studied cylindrically symmet¬ 
ric gravitational collapse with magnetic held and concluded that magnetic 
held can prevent the collapse of cylinder before singularity formation. Ard- 
van and Partovi (1977) investigated dust solution of the held equations with 
electromagnetic held and found that the electrostatic force is balanced by 
gravitational force during collapse of charged dust. The effects of electro¬ 
magnetic held on structure scalars and dynamics of self-gravitating objects 
have been explored by Herrera and his collaborators (Herrera et al. 2011, 
Diprisco et al. 2007). Sharif and his collaborators (Sharif and Bhatti 2012a, 
Sharif and Bhatti 2012b, Sharif and Bhatti 2013a, Sharif and Bhatti 2013b, 
Sharif and Bhatti 2014, Sharif and Yousaf 2012, Sharif and Kausar 2011) 
have extended this work for cylindrical and plane symmetries in GR as well 
as in f(R) gravity with electromagnetic held. 

This paper is organized as follows: In section 2, charged anisotropic cylin¬ 
drical source and Einstcin-Maxwell equations have presented. Section 3 is 
devoted to the generating solutions which represent gravitational collapse 
and expansion of the self-gravitating charged cylinder. We summaries the 
results of the paper in the last section. 

2 Matter Distribution and Field Equations 

This section deals with the interior matter distribution and corresponding 
Einstein-Maxwell’s equations. The non-static spacetime with cylindrical 
symmetry in the interior region of a star is given by (Sharif and Bhatti 
2015) 

ds 2 = — A 2 dt 2 + B 2 dr 2 + C 2 d0 2 + dz 2 , (1) 

where —oo<t<oo, 0 < r < oo, 0 < 0 < 27r, — oo < z < oo are the restric¬ 
tions on the coordinates of cylinder. In the interior region of cylindrically 
symmetric star, we have considered the charged anisotropic fluid for which 
energy-momentum is given by (Di Prisco et al. 2009) 

T a /3 = {/I + Pr)v a Vp ~ ( p r - p z )s a Sp - (p r - p S )k a kp 

+ PrQap + (f2 F P 1 ~ , (2) 

where /j, p r , pe, and p z are the energy density, pressures in r, 9 and z direc¬ 
tions, respectively. Further, v a is four-velocity and s a , k a are four-vectors. 
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Also, F af 3 = — 4*a,f} + 4>p,a is the Maxwell field tensor with four-potential (f> a . 
Moreover, s a and k a are the unit four-vectors which satisfy the following 
relations 


s a s a = k a k a = 1 , 


v a v a = - 1 , 


s Q k n = v a k n = v a s n = 0 . 


In comoving coordinate system, these quantities can be written as 


k a = C5l, v a = 


(3) 


The Maxwell field’s equations are 

F a % = 4tt J a , F [a w = 0, 


where J a is the four-current. In comoving coordinates, the charge inside the 
cylinder is at rest, so we can define the four-potential and four-current as 
follows: 


( fra 


K, 


r 


tv 


where ((r, t) and (f>(r, t) are charge density and scalar potential, respectively. 
The expansion scalar is 


0 


1 B C 
A\B + C 


(4) 


We define the dimensionless anisotropy as follows: 


A a 



(5) 


The corresponding Einstein-Maxwell’s equations have the following form 
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where E = with s(r) = 47T J Q r (BCdr is the total amount of charge per 
unit length of the cylinder. 

Thorne (1965) defined the mass function for cylindrical geometry in the 
form of gravitational C-energy per unit length of the cylinder. The specific 
energy m — El [l is the length of cylinder, i.e., g zz ) of cylindrical geometry 
(fill in the presence of electric charge is given by (Sharif and Bhatti 2015) 


m(t, r) = sC + 


1 
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( 11 ) 


The auxiliary solution of Eq. ((TJ) is 

= B = C\ (12) 

where a is arbitrary constant. Now using Eq. (fT2]) in Eq. (££]), we get the 
following form of expansion scalar 

0 = (l + a)C a ~ l . (13) 


For a > — 1 and a < —1, we obtain expanding and collapsing solutions 
respectively. Using Eq. ffTTIl . in Eqs.(1Sl)- (TTUl) . we get the following form of 
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Einstein-Maxwell’s Equations: 
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For specific values of C(r,t) and a, we can find anisotropic configuration. 
In this case mass function along with electromagnetic field given in Eq. flTTTl 
takes the following form: 


/"V 2 

8m - 8sC - 1 = C 2a - —- 

C 2a 

If C' = C 2a then above equation gives 



(18) 


( 19 ) 


6 


























Figure 1: Density variation for s = 2 and h\{t) — 1 + 1. 


where m > This implies that gravitational collapse leads to the formation 
of a trapping surface at C = ^ (m — |). Also, the integration of trapping 
condition C = C 2a yields 

C^~ 2a ) = (i-2a)r + h(t), (20) 

where h(t) is an arbitrary function. 

3 Generating solution 

For negative and positive values of a , we have collapsing and expanding 
solutions, respectively as follows: 

3.1 Gravitational Collapse with a = — | 

For collapse, expansion scalar will be negative, from Eq. ([TT)h 0 < 0 when 
a > — 1. We assume that a — — | and the condition C' = C 2a , leads to 
C' = C~3, the integration of this equation yields, 

Ctrap — (4r + , (21) 
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Figure 2: Radial pressure variation for s = 2 and hi(t) = 1 + t. 



Figure 3: Transverse pressure variation for s = 2 and h\(t) — 1 + 1. 







Figure 4: Longitudinal pressure variation for s = 2 and h\(t) — 1 + 1. 



Figure 5: Dimensionless anisotropic parameter variation for s = 2 and 
hi(t) — 1 + t. 
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where hi(t) is arbitrary function of time. For a — — |, Eqs. (TT41) - (1T7T) are 
given by 
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Using Eq. 
equations: 
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The dimensionless measure of anisotropy is given by the following equation: 
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For a = 4r, we obtain 0 < 0 and matter density increases for the arbi¬ 
trary choice charge s and time profile hi — 1 + t. As density is increasing 
(see figure. 1), so cylinder goes on collapsing to a point. In this case, the 
length of cylinder is constant, i.e., g zz = 1, therefore base to base collapse is 
impossible, there is only possibility of wall to wall collapse. The anisotropic 
parameter A a changes its sign from negative to positive. The anisotropy will 
be directed outward when pg > p r , this implies that A a > 0 and directed 
inward when pg < p r implying A a < 0. In this case A a > 0, for larger value 
of r as shown in figures 5. This implies that anisotropic force allows the con¬ 
struction of more massive star while A a < 0 near the center, so there exist 
an attractive force. When we talk about an external (electromagnetic) field 
in gravitational, then there is an external force which may distort the generic 
properties of spacetime effectively. So, in the present case electromagnetic 
field enhances the anisotropy and the homogeneity of collapsing star. 
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3.2 Expansion with a = | 


We know that for expansion, the expansion scalar will be positive, from 
Eg. (1111) . 0 < 0, when a > 0. In this case assume that C = (r 2 + rQ) _1 + h 2 (t), 
where /i 2 (f) and are arbitrary function and constant respectively. For 
a — §, Eqs. (fT4]b(H7|) take the following form: 
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If F(t, r) = 1 + h 2 (t)(r 2 + r„) and C = r2 h +r ^ then Eqs. (|32T) - (|34l) become: 
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Figure 6: Density variation for s = 2 and h 2 (t) — 1 + 1. 


The dimensionless measure of anisotropy in this case takes the following form: 


A a = 


3 F 4 
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39) 


For a = §, we have 0 > 0 and matter density decreases (see figure 6) 
for the arbitrary choice of charge and time profile h 2 — 1 + t. In this case p r 
and pe are positive and negative (see figure 7,8), respectively. This implies 
that there exist anisotropy due to opposite behaviour of pressure compo¬ 
nents. From figure 10 A a > 0, there exists repulsive force which causes the 
expansion of matter in this case. The expansion process causes to separate 
the charges apart from each other, this results to weak electromagnetic field 
intensity. 
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Figure 7: Radial pressure variation for s = 2 and h 2 (t) = 1 +t. 



Figure 8: Transverse pressure variation for s = 2 and h 2 {t ) = 1 + t. 
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Figure 9: Longitudinal pressure variation for s — 2 and h 2 {t) — 1 + 1. 



Figure 10: Dimensionless anisotropic parameter variation for s = 2 and 
h 2 (t) — 1 + t. 
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4 Concluding Remarks 

During the last few decades, there has been a growing interest to study the 
relativistic anisotropic systems due to the existence of such systems in as¬ 
tronomical objects. The exact solutions of anisotropic sources are helpful 
to determine the anisotropy of the universe during any era. The effects of 
anisotropy on the late-time expansion of inhomogeneous universe have been 
studied by Barrow and Maartens (1998). They remarked that the decrease 
in shear anisotropy can be figure out by measuring the anisotropic pressure 
of the cosmological model. Herrera and Santos (1997) pointed out that the 
phase transition in anisotropic highly dense system would occur during the 
gravitational collapse. Further, they concluded that such system may tran¬ 
sited to a pion condensed phase, where a soften equation of state can provide 
enough exhausted energy. 

This paper is aimed to study the generating solution of Einstein-Maxwell 
held equations with anisotropic cylindrically symmetric fluid. We have used 
the auxiliary solution of one held equation to determine the solution of the re¬ 
maining equations. The application of assumed solution in expansion scalar, 
allows us to determine the range of free constant a, for which expansion scalar 
0 is positive or negative, leading to expansion and collapse. The C-energy 
analogous to Misner-Sharp mass has been calculated with the contribution. 
We have impose the condition C' = C 2a , on the mass function which leads 
to the existence of trapping horizon at C — ^ {m — |), provided m > |. In 
this case curvature singularity is hidden at the center of trapping horizon. 

The expansion scalar 0 = ( 2+a)C^ a \ becomes 0 = 0, for a = —1, 0 > 
0, for a > —1,0 < 0, for a < —1, which corresponds to bouncing, expansion 
and collapse, respectively. In other words 0 > 0 , a G (0, oo), and 0 < 0 for 
a G (—oo, 0) which corresponds to expansion and collapse, respectively. For 
the sake of simplicity, we have taken a = — § for gravitational collapse and 
a — | for expansion, explicitly. The full dynamics of the system has been 
discussed in both cases. The matter density is increasing/dcreasing function 
with arbitrary choice of charge parameter and time profiles. The pressures 
p r . po and p z are different in both cases, therefore the pressure anisotropy is 
non-vanishing in both cases. This anisotropy is increasing function in both 
cases, it is due to presence of electromagnetic held, as it produces an external 
repulsive force to distort the geometry of the star. 

We would like to mention this work with charged plane symmetric 
source is in progress. 
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